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Q\ • Abstract 

Recent progress in solar Hard X-ray (HXR) observations with RHESSI data and 
^ ■ methods for spectral inversion allow us to study model-independent mean electron 

OO , flux spectra in solar flares. We report several hard X-ray events observed by RHESSI 

in which the photon spectra /(e) are such that the inferred source mean electron 
00 ' spectra are not consistent with the standard model of collisional transport in solar 

, flares. The observed photon spectra are so flat locally that the recovered mean 

(<~^ I electron flux spectra show a dip around 17-31 keV. While we note that alternative 

' explanations, unrelated to electron transport, have not been ruled out, we focus 

Q^. on the physical implications of this tentative result for the collisional thick-target 

Q ! model. 
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1 Introduction 



It has long been recognized (Brown, 1971; Brown and Emslie, 1988; Johns 
and Lin, 1992) that even spatially integrated spectra /(e) of flare hard X-ray 
(HXR) collisional bremsstrahlung bursts carry crucial information on flare 
electron acceleration, transport, and energy budget. Even the absence of spec- 
tral features, such as in an (energy-scale-free) pure power-law, is likely to in- 
dicate stochastic multi-scale processes. On the other hand, any features which 
exist in the acceleration spectrum are smeared out in the propagation and 
radiation process so that any features detected in /(e) have potentially very 
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strong implications for the propagation (Kontar ct al, 2003) and accelera- 
tion processes, and for the electron energy budget in terms of any low energy 
cut-off. 



The large data base of high resolution spectra from the RHESSI mission (Lin 
et al, 2002) has made possible a thorough search for such special HXR spec- 
tral diagnostic properties. Care has to be taken, however, not to misinterpret 
features as entirely real solar ones, such as the 50 keV electron spectral feature 
in the July 23, 2002 event (Plana et al, 2003). This can be attributed partly 
to pulse pile- up at the high count rate in this event. Two approaches are used 
in interpreting 7(e) data. One is to adopt a parametric form for the mean 
radiation source electron spectrum (or for the electron 'injection' spectrum ) 
and predict /(e), using an accurate bremsstrahlung cross-section Q{e, E), and 
search for the best fit in the model parameter space. This is commonly done 
using a source-mean flux spectra F{E) containing an isothermal Maxwellian, 
Ft{E), defined by an emission measure EM and temperature T, plus a non- 
thermal Fpow{E) parameterized as a single or double power-law with a sharp 
low energy cut-off, Ec (Holman et al, 2003). Such best fits tend to yield an 
i?c ~ 10 — 30 keV, well above kT ~ 1 — 3 keV, but with EM large enough 
that Ec is 'buried' in the Ft{E) 'tail' and it is impossible to be sure whether 
Fpow{E) has an actual cut off or blends smoothly with Ft{E), especially if 
there is some spread in T. Recently, Schwartz et al (2003) have reported, for 
the August 20, 2002 fiare that a best fit requires an E^. 30 keV in FpowiE) 
seen clearly above the isothermal best fit Fx{E). To evaluate whether such 
a real feature/cut-off is demanded by the data, as opposed to simply being 
compatible with it, it is essential to adopt a non-parametric approach to in- 
terpreting /(e) - i.e. to infer from /(e) what range of functions F{E) allows a 
statistically acceptable fit to /(e). This inverse/inferential approach has been 
adopted, developed and applied to data by a number of authors (Craig and 
Brown 1986, Johns and Lin, 1992; Thomson et al., 1992; Plana et al, 2003; 
Kontar et al, 2004; Kontar et al, 2005). 



In this paper we apply one of the best available inversion algorithms of Kontar 
et al (2004) to several events we have found in the RHESSI database to show 
real gaps or dips in F{E) probably demanded by, and not just consistent 
with, the data. The proposed implications of these results for flare electron 
acceleration, propagation, and energy budget are briefly discussed. 
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2 Bremsstrahlung Source Models and Constraints on the Emission 
Spectra 



2.1 Mean source electron spectra F{E) 

In this section we generalize to arbitrary cross-section Q{e,E) discussions of 
bremsstrahlung spectrum constraints for models discussed earlier by various 
authors. 

Following Brown (1971) and Brown, Emshe and Kontar (2003) we emphasize 
that the electron distribution function which can be inferred from 7(e) without 
source model assumptions (apart from optical thinness and isotropy) is the 
density weighted mean radiation source electron flux spectrum F{E) (electrons 
cm~^ s~^ keV~^) defined by 



where F{E,r), n(r) are the local electron flux spectrum and source proton 
density at position r in radiating volume V with mean target proton density 
n = J n(r) dV. F{E) uniquely defines the bremsstrahlung spectrum at 
the earth by 



where Q{e,E) is the cross-section differential in e (Haug, 1997). In general 

Q{e, E) should include angle averaging to allow for anisotropy of the electron 
distribution function. The only correction we do not discuss explicitly here is 
for photospheric back scatter, though it might be important for our results. 

The general approach to Equation (2) is to treat it as an integral equation to 
reconstruct F{E) from /(e) by deconvolution through Q{e, E) using regulariza- 
tion techniques e.g. (Craig and Brown, 1986; Kontar et al, 2004). Limitations 
can be placed on 7(e) for a physically acceptable solution of (2) to exist. In 
particular, generalizing (Brown and Emshe, 1988), 
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Since any physical F{E) should be nonnegativc and since d{Q{e, E)e)/de < 
Ve > for the relativistic cross-section of Haug, (1997), we immediately have 
that d{I{e)e)/de < 0, Ve > 0. It is useful to consider the shape of /(e) in terms 
of the local spectral index 7(e) defined by 

7(e) = (4) 
lie) de ' ^ ^ 



(note that definition (4) does not correspond to /(e) ~ e^'''*^'^) except for con- 
stant 7 (Conway et al, 2003)). Then in terms of the spectral index the condition 
d{I{e)e)/de < can be expressed 

7(e) > 1, Ve (5) 



i.e. any bremsstrahlung emission spectrum /(e) should decrease with a loga- 
rithmic gradient larger than 1. If the condition (5) is violated then the spec- 
trum /(e) is not from an optically thin bremsstrahlung source. 

To date no HXR spectrum from a solar fiare violating Condition (5) has ever 
been reported, so there is no challenge to the belief that fiare HXR burst spec- 
tra are consistent with optically thin bremsstrahlung. However, more stringent 
spectral compatibility conditions apply to specific models of how F{E) is 
formed from an injection spectrum Fq{Eq). 



2.2 Thick target model 



A thick target is a source in which electrons injected into the source decelerate 
to rest under the action of energy loss processes, within the observed integra- 
tion time and volume. The terminology 'thick' is similar to optical depth, a 
thick target being a medium where plasma electron column density is suffi- 
cient to stop energetic electrons, or, equally, the mean free path of the fast 
electrons is much less than the system length. In other words, for energy loss 
cross section Qe = —dlog{E)/dN, NQe 3> 1 within the column density 
of V at all E of interest. The resulting mean electron fiux spectrum is then 
the density-weighted average of the fiux spectrum of the electrons along their 
energy loss path z, as governed by Qe for all the energy loss processes involved. 

The assumptions in inferring Fq{Eq) from the total /(e) are that : (a) the 
bremsstrahlung emission from acceleration region can be neglected, and (b) 
physically distinct 'acceleration' and 'target' (deceleration) regions exist. 

Let us consider the motion for an electron in the target from an injection 
spectrum Fq{Eq) {E{z = 0) = £"0) for a given background plasma density 
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profile n{z), and introdTicc the cohimn depth N{z) of the plasma so that 
n{z) = dN/dz. It is important to note that z is an electron path. Therefore, a 
column depth can be also defined as dN{z) = n{z)v{z)dt. Comparing with the 
original thick target model (Brown, 1971) where only coUisional losses have 
been taken into account we get, for that case 

dE , , , 

— = -K E 6 
dNcoii ' ^ ' 



where K — 27re^A, e is the electron charge,and A f=:i 20 — 25 is the Coulomb 
logarithm. 

Assuming that we can find the solution of the equation of motion of an electron 
in the form E — E(Eq, N) we can write down the mean electron flux in the 
following form 

A °° 



F{E)^— I F{E{Eo,N))dN. (7) 







where A is the injected area. The evolving electron spectrum F{E{Eo, N)) 
can then be expressed in terms of the initial or injected electron spectrum 
using the continuity equation F{E{Eq, N))dE = Fo{EQ)dEo. After the change 
of variables N ^ Eo in (7) one has for the mean electron spectrum with 
arbitrary energy losses. 

m--^ J MEo)[^) dE.. (8) 

E ^ ^ 



(The mean electron flux given by (8) makes use of deterministic electron prop- 
agation ignoring collective and diffusion effects). In previous papers (Brown 
and MacKinnon, 1985) it was also assumed that the upper limit in the integral 
(8) is infinity ^ , e.g. Eo{N oo, E) = oo, which is true for purely coUisional 
losses. Equations (6) and (8) become 

A TP c 

F{E) = J Fo{Eo)dEo. (9) 

E 



^ Actually Eo{N ^ oo, £^) = oo for any non-collisional losses described by the 
equation of the type ^ = —(p{E) 
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2.3 Purely collisional thick target 



We now consider the case of a purely collisional thick target (9), when the 
injection spectrum takes the form 



KnV d F{E) 



A dE E 

.nV F 



E=Eo 



K^^{i+m)E=Eo (10) 



where S{E) — —dlnF/dlnE is the spectral index of the density weighted 
mean electron flux distribution. Prom equation (10) it follows that, for a pho- 
ton spectrum /(e) to be produced by a purely collisional thick target with a 
physical (nonnegative) injection spectrum Fq{Eq) > 0, the mean electron fiux 
F{E) derived from the photon spectrum /(e) must have a logarithmic slope 

5(E) > -1 V£; ^ (11) 



(or dlogF/dlogE < +1). 

In the next section we report RHESSl /(e) data in which the above condition 
appears to be violated, and a purely collisional thick target interpretation is 
impossible. 

Equation (9) also explicitly says that if the injected spectrum Fq{Eq) has a 
low energy cut-off at energy Eq, then the purely collisional thick target mean 
spectrum F{E) can have a spectral index S{E) — —1, for £■ < Eq i.e. a 
logarithmic slope of > -|-1) which is the lowest S value that can be explained 
by a collisional thick target. However, if the spectral index at some point is 
less than -1, then condition (11) is violated, and non-coUisional models must 
be involved. 



3 RHESSl observations 



There are flares observed by RHESSl (Lin et al, 2002) with unusually fiat pho- 
ton spectra resulting in peculiar mean electron distribution functions F{E). 
Unhke the July 23, 2002 flare (Plana et al, 2003), in these flares the count rate 
in the front segments of the RHESSl detectors were too low for pulse pile-up 
(Smith et al, 2002) to be important. Forward fits applied to this data show 
that the X-ray spectra can be best fitted by a broken power-law with a low 
energy cutoff in the range between 15 and 35 keV. This suggests that there 
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are no injected nonthermal electrons in the range below a few tens of kcV. 
To verify this we used a recently developed rcgularization algorithm (Kontar 
et al. 2004) to infer model - free electron spectra for these events. The mean 
electron flux spectrum inferred shows a clear dip at energies between 15 and 
35 keV. A similar feature, though less clear and at the energies around 50keV, 
has been recovered by Plana et al (2003). The main results are presented in 
Fig 1-3. 

3.1 Inversion of photon spectra 

The observed spectrum of the flare is the convolution of the cross-section and 
mean electron flux. The problem of inferring the mean electron spectrum is ill- 
posed (Craig and Brown, 1986). Thus, one has to avoid unphysical behaviour 
in F[E) using some constraints. This process is called regularization. 

To infer the mean electron spectrum we used the following recently developed 
regularization algorithm of Kontar et al, (2004) 

C(¥) = II AF - If + A||LFf , £(F) = min (12) 

where A is a regularization parameter, L is a regularizing first order derivative 
operator, and A is the matrix representation of the cross-section operator A 
in 

1 °° 

{AF){e) = ^nvl F{E) Q{e, E) dE, (13) 

e 

i.e. comprises binned integral components J^^^ Q{^i, E)dE. For our calcula- 
tions we used the Haug (1997) cross-section. 

One shonid note that if wc try to infer the electron spectra without additional 
constraints, e.g. setting to zero the penalty term in Eq. 12, the resulting F{E) 
would be a violently oscillating function. Such unphysical results come from 
the ill-posedness of the problem, since small noise perturbations in /(e) will 
be strongly amplified. 

3.2 Mean Electron Spectrum and Spectral Index 

For our analysis we have selected three flares with unusually flat local spectra. 
The data reduction has been performed in the same way as in Kontar et al 
(2005). Then the regularization algorithm has been applied to flnd the mean 
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Fig. 1. Spatially integrated photon flux /(e) (upper panel), mean electron flux F{E) 
and electron spectral index d{E) (lower panel) for 25 April 2002 05:55 UT flare 
using first order regularization. The photon spectrum has been accumulated over 
the impulsive phase of the flare. 



electron flux and its spectral index. Figures 1-3 show fiat mean electron flux 
spectra with a positive derivative in the range 15-35 keV. The depth of the 
dip and the resulting electron spectral index is different from flare to flare, 
though qualitatively similar. The April 25 event (Figure 1) has a dip at 19 
keV but the error bars do not allow us to conclude that this dip is real. 

The confidence intervals on the regularized solution and spectral index (Fig- 
ures 1-3) have been calculated as a maximum deviation from the solution in 
a set of 30 random realizations of photon data within a ilcr range. The dis- 
tribution of errors has been discussed in Plana et al 2003; Kontar et al, 2004. 
To allow for instrumental uncertainties, la was taken to be not less than 3% 
of the photon flux, and thus can be treated as an upper estimate of the error. 

The August 2002 (Figure 2) flare shows an extremely clear gap around 20 
keV in the mean electron flux. Forward-fit to the data shows an unambiguous 
low energy cut-off at ~ 31 keV (Kasparova et al, 2005). To verify the reality 
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Fig. 2. The same as in Fig 1 for 20 August 2002 08:24 UT flare. 

of this feature, we used different orders of regularization for the inversion 
of the corresponding flare X-ray spectra. All zero, first, and second order 
regularization methods produce the dip (Figure 2). 

Figure 3 also shows that at least three successive points have S less than — 1 
with la confidence and 4 points less than zero. The former suggests that the 
thick-target model is unacceptable while the latter indicates a low energy cut- 
off in the injected electron spectrum. Though we have used the la criterion, 
the result is much more significant than that 65% of a single la deviation since 
there are three successive such deviations. While a more rigorous statistical 
analysis is necessary, this should be significant at roughly the level of 1 — (1 — 
0.65)^ or about 97% which makes it highly suggestive at least. 



4 Discussion and conclusions 



Our analysis of the observed photon spectra shows features that have not been 
observed before. Johns and Lin (1992) reported the downturn of the derived 
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Fig. 3. The same as in Fig 1 for 17 Sept 2002 05:55 UT flare. 

electron spectrum for earlier balloon data, though the errors were too big to be 
conclusive. In the results presented in this paper we are much more confident 
that a dip in the electron distribution is required by the photon spectra, though 
we caution that the photon spectrum itself is subject to additional corrections 
for photospheric albedo, directivity, and instrument effects (Smith et al, 2002). 



One of the important implications of our result if upheld is related to the accel- 
erated spectrum Fq^Eq). As can be seen from the results, the photon spectra 
require a low energy cut-off in the injected spectra Fq^Eq) to be produced, e.g. 
the absence of electrons in the range below 20-30 keV. This is quite a strong 
requirement on acceleration models, that often produce extended power law 
spectra above thermal energies of a few keV. The cut-off poses physical chal- 
lenges to electron acceleration and transport since an electron beam with a 
velocity distribution function such that f{v) ~ F{v'^) should be unstable to 
the generation of plasma waves if df{v)/dv > 0. 

We stress that, if the existence of a low energy cut-off in the injected electron 
spectrum Fq^Eq) is confirmed, its implications for total fiare electron power, 
acceleration and propagation are profound. Moreover, the evidence of a steep 
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turn over in F{E) would reject the long-standing collisional transport thick- 
target model of flare electron propagation (Brown, 1971) where the accelerated 
electron spectrum Fq{Eq) is modified solely by Coulomb collisions along the 
electron path. 
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